Search for: All records

Observations indicate dust populations vary between galaxies and within them, suggesting a complex life cycle and evolutionary history. Here we investigate the evolution of galactic dust populations across cosmic time using a suite of cosmological zoom-in simulations from the Feedback in Realistic Environments project, spanning $M_{\rm vir}=10^{9-12}{M}_{\odot };\, M_{*}=10^{6-11}\, {M}_{\odot }$. Our simulations incorporate a dust evolution model that accounts for the dominant sources of dust production, growth, and destruction and follows the evolution of specific dust species. All galactic dust populations in our suite exhibit similar evolutionary histories, with gas–dust accretion being the dominant producer of dust mass for all but the most metal-poor galaxies. Similar to previous works, we find the onset of efficient gas–dust accretion occurs above a ‘critical’ metallicity threshold (Zcrit). Due to this threshold, our simulations reproduce observed trends between galactic D/Z and metallicity and element depletion trends in the interstellar medium. However, we find Zcrit varies between dust species due to differences in key element abundances, dust physical properties, and life cycle processes resulting in $Z_{\rm crit}\sim 0.05{\rm Z}_{\odot },\, 0.2{\rm Z}_{\odot },\, 0.5{\rm Z}_{\odot }$ for metallic iron, silicates, and carbonaceous dust, respectively. These variations could explain the lack of small carbonaceous grains observed in the Magellanic Clouds. We also find a delay between the onset of gas–dust accretion and when a dust population reaches equilibrium, which we call the equilibrium time-scale (τequil). The relation between τequil and the metal enrichment time-scale of a galaxy, determined by its recent evolutionary history, can contribute to the scatter in the observed relation between galactic D/Z and metallicity.

 

We measure the CO-to-H₂conversion factor (α_CO) in 37 galaxies at 2 kpc resolution, using the dust surface density inferred from far-infrared emission as a tracer of the gas surface density and assuming a constant dust-to-metal ratio. In total, we have ∼790 and ∼610 independent measurements ofα_COfor CO (2–1) and (1–0), respectively. The mean values forα_{CO (2–1)}andα_{CO (1–0)}are${9.3}_{-5.4}^{+4.6}$and${4.2}_{-2.0}^{+1.9}{M}_{\odot }{\mathrm{pc}}^{-2}{(\mathrm{K}\mathrm{km}{\mathrm{s}}^{-1})}^{-1}$, respectively. The CO-intensity-weighted mean is 5.69 forα_{CO (2–1)}and 3.33 forα_{CO (1–0)}. We examine howα_COscales with several physical quantities, e.g., the star formation rate (SFR), stellar mass, and dust-mass-weighted average interstellar radiation field strength ($\overline{U}$). Among them,$\overline{U}$, Σ_SFR, and the integrated CO intensity (W_CO) have the strongest anticorrelation with spatially resolvedα_CO. We provide linear regression results toα_COfor all quantities tested. At galaxy-integrated scales, we observe significant correlations betweenα_COandW_CO, metallicity,$\overline{U}$, and Σ_SFR. We also find thatα_COin each galaxy decreases with the stellar mass surface density (Σ_⋆) in high-surface-density regions (Σ_⋆≥ 100M_⊙pc⁻²), following the power-law relations${\alpha }_{\mathrm{CO}(2–1)}\propto {\mathrm{\Sigma }}_{\star }^{-0.5}$and${\alpha }_{\mathrm{CO}(1–0)}\propto {\mathrm{\Sigma }}_{\star }^{-0.2}$. The power-law index is insensitive to the assumed dust-to-metal ratio. We interpret the decrease inα_COwith increasing Σ_⋆as a result of higher velocity dispersion compared to isolated, self-gravitating clouds due to the additional gravitational force from stellar sources, which leads to the reduction inα_CO. The decrease inα_COat high Σ_⋆is important for accurately assessing molecular gas content and star formation efficiency in the centers of galaxies, which bridge “Milky Way–like” to “starburst-like” conversion factors.

 

Young stellar objects (YSOs) are the gold standard for tracing star formation in galaxies but have been unobservable beyond the Milky Way and Magellanic Clouds. But that all changed when the JWST was launched, which we use to identify YSOs in the Local Group galaxy M33, marking the first time that individual YSOs have been identified at these large distances. We present Mid-Infrared Instrument (MIRI) imaging mosaics at 5.6 and 21 $\mu$m that cover a significant portion of one of M33’s spiral arms that has existing panchromatic imaging from the Hubble Space Telescope and deep Atacama Large Millimeter/submillimeter Array CO measurements. Using these MIRI and Hubble Space Telescope images, we identify point sources using the new dolphot MIRI module. We identify 793 candidate YSOs from cuts based on colour, proximity to giant molecular clouds (GMCs), and visual inspection. Similar to Milky Way GMCs, we find that higher mass GMCs contain more YSOs and YSO emission, which further show YSOs identify star formation better than most tracers that cannot capture this relationship at cloud scales. We find evidence of enhanced star formation efficiency in the southern spiral arm by comparing the YSOs to the molecular gas mass.

 

We present a method to characterize star-formation driven outflows from edge-on galaxies and apply this method to the metal-poor starburst galaxy, Mrk 1486. Our method uses the distribution of emission line flux (from H β and [O iii] 5007) to identify the location of the outflow and measure the extent above the disc, the opening angle, and the transverse kinematics. We show that this simple technique recovers a similar distribution of the outflow without requiring complex modelling of line-splitting or multi-Gaussian components, and is therefore applicable to lower spectral resolution data. In Mrk 1486 we observe an asymmetric outflow in both the location of the peak flux and total flux from each lobe. We estimate an opening angle of 17–37° depending on the method and assumptions adopted. Within the minor axis outflows, we estimate a total mass outflow rate of ∼2.5 M⊙ yr−1, which corresponds to a mass loading factor of η = 0.7. We observe a non-negligible amount of flux from ionized gas outflowing along the edge of the disc (perpendicular to the biconical components), with a mass outflow rate ∼0.9 M⊙ yr−1. Our results are intended to demonstrate a method that can be applied to high-throughput low spectral resolution observations, such as narrow-band filters or low spectral resolution integral field spectrographs that may be more able to recover the faint emission from outflows.

 

Carbon monoxide (CO) emission constitutes the most widely used tracer of the bulk molecular gas in the interstellar medium (ISM) in extragalactic studies. The CO-to-H 2 conversion factor, α 12 CO(1−0) , links the observed CO emission to the total molecular gas mass. However, no single prescription perfectly describes the variation of α 12 CO(1−0) across all environments within and across galaxies as a function of metallicity, molecular gas opacity, line excitation, and other factors. Using spectral line observations of CO and its isotopologues mapped across a nearby galaxy, we can constrain the molecular gas conditions and link them to a variation in α 12 CO(1−0) . Here, we present new, wide-field (10 × 10 arcmin 2 ) IRAM 30-m telescope 1 mm and 3 mm line observations of 12 CO, 13 CO, and C 18 O across the nearby, grand-design, spiral galaxy M101. From the CO isotopologue line ratio analysis alone, we find that selective nucleosynthesis and changes in the opacity are the main drivers of the variation in the line emission across the galaxy. In a further analysis step, we estimated α 12 CO(1−0) using different approaches, including (i) via the dust mass surface density derived from far-IR emission as an independent tracer of the total gas surface density and (ii) local thermal equilibrium (LTE) based measurements using the optically thin 13 CO(1–0) intensity. We find an average value of ⟨ α 12 CO(1 − 0) ⟩ = 4.4  ±  0.9  M ⊙  pc −2  (K km s −1 ) −1 across the disk of the galaxy, with a decrease by a factor of 10 toward the 2 kpc central region. In contrast, we find LTE-based α 12 CO(1−0) values are lower by a factor of 2–3 across the disk relative to the dust-based result. Accounting for α 12 CO(1−0) variations, we found significantly reduced molecular gas depletion time by a factor 10 in the galaxy’s center. In conclusion, our result suggests implications for commonly derived scaling relations, such as an underestimation of the slope of the Kennicutt Schmidt law, if α 12 CO(1−0) variations are not accounted for. 

Determining how the galactic environment, especially the high gas densities and complex dynamics in bar-fed galaxy centers, alters the star formation efficiency (SFE) of molecular gas is critical to understanding galaxy evolution. However, these same physical or dynamical effects also alter the emissivity properties of CO, leading to variations in the CO-to-H₂conversion factor (α_CO) that impact the assessment of the gas column densities and thus of the SFE. To address such issues, we investigate the dependence ofα_COon the local CO velocity dispersion at 150 pc scales using a new set of dust-basedα_COmeasurements and propose a newα_COprescription that accounts for CO emissivity variations across galaxies. Based on this prescription, we estimate the SFE in a sample of 65 galaxies from the PHANGS–Atacama Large Millimeter/submillimeter Array survey. We find increasing SFE toward high-surface-density regions like galaxy centers, while using a constant or metallicity-basedα_COresults in a more homogeneous SFE throughout the centers and disks. Our prescription further reveals a mean molecular gas depletion time of 700 Myr in the centers of barred galaxies, which is overall three to four times shorter than in nonbarred galaxy centers or the disks. Across the galaxy disks, the depletion time is consistently around 2–3 Gyr, regardless of the choice ofα_COprescription. All together, our results suggest that the high level of star formation activity in barred centers is not simply due to an increased amount of molecular gas, but also to an enhanced SFE compared to nonbarred centers or disk regions.

 

Abstract NGC 602 is a young, low-metallicity star cluster in the “Wing” of the Small Magellanic Cloud. We reveal the recent evolutionary past of the cluster through analysis of high-resolution (∼0.4 pc) Atacama Large Millimeter/submillimeter Array observations of molecular gas in the associated H ii region N90. We identify 110 molecular clumps ( R < 0.8 pc) traced by CO emission, and study the relationship between the clumps and associated young stellar objects (YSOs) and pre-main-sequence (PMS) stars. The clumps have high virial parameters (typical α vir = 4–11) and may retain signatures of a collision in the last ≲8 Myr between H i components of the adjacent supergiant shell SMC-SGS 1. We obtain a CO-bright-to-H 2 gas conversion factor of X CO, B = (3.4 ± 0.2) × 10 20 cm −2 (K km s −1 ) −1 , and correct observed clump properties for CO-dark H 2 gas to derive a total molecular gas mass in N90 of 16,600 ± 2400 M ⊙ . We derive a recent (≲1 Myr) star formation rate of 130 ± 30 M ⊙ Myr −1 with an efficiency of 8% ± 3% assessed through comparing total YSO mass to total molecular gas mass. Very few significant radial trends exist between clump properties or PMS star ages and distance from NGC 602. We do not find evidence for a triggered star formation scenario among the youngest (≲2 Myr) stellar generations, and instead conclude that a sequential star formation process in which NGC 602 did not directly cause recent star formation in the region is likely. 

Abstract The CO-to-H 2 conversion factor ( α CO ) is central to measuring the amount and properties of molecular gas. It is known to vary with environmental conditions, and previous studies have revealed lower α CO in the centers of some barred galaxies on kiloparsec scales. To unveil the physical drivers of such variations, we obtained Atacama Large Millimeter/submillimeter Array bands (3), (6), and (7) observations toward the inner ∼2 kpc of NGC 3627 and NGC 4321 tracing 12 CO, 13 CO, and C 18 O lines on ∼100 pc scales. Our multiline modeling and Bayesian likelihood analysis of these data sets reveal variations of molecular gas density, temperature, optical depth, and velocity dispersion, which are among the key drivers of α CO . The central 300 pc nuclei in both galaxies show strong enhancement of temperature T k ≳ 100 K and density n H 2 > 10 3 cm −3 . Assuming a CO-to-H 2 abundance of 3 × 10 −4 , we derive 4–15 times lower α CO than the Galactic value across our maps, which agrees well with previous kiloparsec-scale measurements. Combining the results with our previous work on NGC 3351, we find a strong correlation of α CO with low- J 12 CO optical depths ( τ CO ), as well as an anticorrelation with T k . The τ CO correlation explains most of the α CO variation in the three galaxy centers, whereas changes in T k influence α CO to second order. Overall, the observed line width and 12 CO/ 13 CO 2–1 line ratio correlate with τ CO variation in these centers, and thus they are useful observational indicators for α CO variation. We also test current simulation-based α CO prescriptions and find a systematic overprediction, which likely originates from the mismatch of gas conditions between our data and the simulations. 

Recent strides have been made developing dust evolution models for galaxy formation simulations but these approaches vary in their assumptions and degree of complexity. Here, we introduce and compare two separate dust evolution models (labelled ‘Elemental’ and ‘Species’), based on recent approaches, incorporated into the gizmo code and coupled with fire-2 stellar feedback and interstellar medium physics. Both models account for turbulent dust diffusion, stellar production of dust, dust growth via gas-dust accretion, and dust destruction from time-resolved supernovae, thermal sputtering in hot gas, and astration. The ‘Elemental’ model tracks the evolution of generalized dust species and utilizes a simple, ‘tunable’ dust growth routine, while the ‘Species’ model tracks the evolution of specific dust species with set chemical compositions and incorporates a physically motivated, two-phase dust growth routine. We test and compare these models in an idealized Milky Way-mass galaxy and find that while both produce reasonable galaxy-integrated dust-to-metals (D/Z) ratios and predict gas-dust accretion as the main dust growth mechanism, a chemically motivated model is needed to reproduce the observed scaling relation between individual element depletions and D/Z with column density and local gas density. We also find the inclusion of theoretical metallic iron and O-bearing dust species are needed in the case of specific dust species in order to match observations of O and Fe depletions, and the integration of a sub-resolution dense molecular gas/CO scheme is needed to both match observed C depletions and ensure carbonaceous dust is not overproduced in dense environments.

 

Abstract The metallicity and gas density dependence of interstellar depletions, the dust-to-gas (D/G), and dust-to-metal (D/M) ratios have important implications for how accurately we can trace the chemical enrichment of the universe, either by using FIR dust emission as a tracer of the ISM or by using spectroscopy of damped Ly α systems to measure chemical abundances over a wide range of redshifts. We collect and compare large samples of depletion measurements in the Milky Way (MW), Large Magellanic Cloud (LMC) ( Z = 0.5 Z ⊙ ), and Small Magellanic Cloud (SMC) ( Z = 0.2 Z ⊙ ). The relations between the depletions of different elements do not strongly vary between the three galaxies, implying that abundance ratios should trace depletions accurately down to 20% solar metallicity. From the depletions, we derive D/G and D/M. The D/G increases with density, consistent with the more efficient accretion of gas-phase metals onto dust grains in the denser ISM. For log N (H) > 21 cm −2 , the depletion of metallicity tracers (S, Zn) exceeds −0.5 dex, even at 20% solar metallicity. The gas fraction of metals increases from the MW to the LMC (factor 3) and SMC (factor 6), compensating for the reduction in total heavy element abundances and resulting in those three galaxies having the same neutral gas-phase metallicities. The D/G derived from depletions are respective factors of 2 (LMC) and 5 (SMC) higher than the D/G derived from FIR, 21 cm, and CO emission, likely due to the combined uncertainties on the dust FIR opacity and on the depletion of carbon and oxygen.